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ABSTRACT 

We examine far-infrared and submillimeter spectral energy distributions for galaxies in the Infrared Space 
Observatory Atlas of Bright Spiral Galaxies. For the 71 galaxies where we had complete 60 - 180 fim data, 
we fit blackbodies with A -1 emissivities and average temperatures of 31 K or A -2 emissivities and average 
temperatures of 22 K. Except for high temperatures determined in some early-type galaxies, the temperatures 
show no dependence on any galaxy characteristic. For the 60 - 850 fim range in eight galaxies, we fit blackbodies 
with A -1 , A~ 2 , and A _/3 (with (3 variable) emissivities to the data. The best results were with the A - ' 3 
emissivities, where the temperatures were ~ 30 K and the emissivity coefficient (3 ranged from 0.9 to 1.9. 
These results produced gas to dust ratios that ranged from 150 to 580, which were consistent with the ratio 
for the Milky Way and which exhibited relatively little dispersion compared to fits with fixed emissivities. 
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1. Introduction 

The far-infrared emission from galaxies is emitted by 
interstellar dust heated by a variety of energy sources. 
Two important energy sources are the general interstel- 
lar radiation field and the more intense radiation in star 
formation regions. These should produce dust sources 
at different dust temperatures and with different spa- 
tial distributions. Understanding these dust tempera- 
ture components is necessary for understanding the rela- 
tion between far-infrared emission and star formation in 
galaxies. 

While much of the research in this subject has focused 
on ultraluminous galaxies and high redshift objects, some 
studies have investigated dust temperatures in normal 
spiral galaxies. However, no clear consensus has been 
reached on the temperatures of the dust components in 
spiral galaxies. Perhaps the most contentious issue is 
whether a significant dust component at a very low tem- 
perature ~ 10-15 K exists. This issue needs to be resolved 
to understand how the interstellar radiation field and star 
formation regions heat dust to produce far-infrared emis- 
sion. 

Before the launch of the Infrared Space Observatory 
(ISO) (Kessler et al. 1996), most research on dust tem- 
peratures in galaxies relied on combining IRAS or Kuipcr 
Airborne Observatory (KAO) data with ground-based 
submillimeter or millimeter data. The launch of ISO 
has permitted better far-infrared spectral energy distri- 
butions to be measured. 
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Smith (1982) was the first to measure a dust temper- 
ature in a non-active spiral. Using KAO data, he found 
a temperature of 20 K in NGC 5194 (M 51). Chini et 
al. (1986) combined submillimctcr bolometer data from 
the NASA IRTF with IRAS far-infrared data to conclude 
that spiral galaxies have cold dust components with an 
average temperature of 16 K. Stark et al. (1989) mea- 
sured far-infrared fluxes with the KAO and submillimeter 
fluxes with the NASA IRTF for three Virgo cluster spi- 
rals and inferred significantly warmer dust temperatures 
than those determined by Chini et al. (1986). Eales, 
Wynn- Williams, & Duncan (1989) measured dust tem- 
peratures in 11 galaxies using IRAS far-infrared data and 
submillimctcr bolometer data from UKIRT and deter- 
mined that the coolest dust component in spiral galaxies 
was 30 - 50 K, much warmer than previous measure- 
ments. These studies thus established a controversy over 
the apparently simple issue of the dust temperatures in 
spiral galaxies. 

Subsequent studies have produced similarly conflict- 
ing results. Most have focused on individual galaxies, 
including NGC 660, UGC 3490 (Chini & Kriigel 1993), 
NGC 3267 (Sievers et al. 1994), NGC 4631 (Braine et al. 
1995), M 51 (Guelin et al. 1995; Hippelein et al. 1996), 
M 101 (Hippelein et al. 1996), the Andromeda Galaxy 
(Haas et al. 1998), and NGC 891 (Popescu et al. 2000). 
These studies typically find dust temperatures of ~ 20 K, 
with the lowest temperature being 16 K and the highest 
being 33 K. There have been a few studies of larger sam- 
ples. Some of these surveys have found far-infrared dust 
temperatures as high as 35 K (Clements, Andrcani, & 
Chase 1993; Dunne et al. 2000). Others have found tem- 
peratures between 20 and 30 K Dunne & Eales (2001); 
Popescu et al. (2002), and still others have found dust 
components as cold as 10-15 K (Chini et al. 1995; Kriigel 
et al. 1998; Siebenmorgen, Kriigel, & Chini 1999; Stickel 
et al. 2000). 

Part of the reason no consensus has been reached on 
dust temperatures in normal galaxies is because of poor 
sampling of the spectral energy distribution (SED). To 
find reliable dust temperatures it is vital to evenly sam- 
ple the SEDs from the Wien side to the Rayleigh-Jeans 
side of the Planck spectrum. Some studies have had such 
poor spectral sampling that Planck functions are fit to 
data with no degrees of freedom, and some have entire 
dust temperature components defined by a single mea- 
surement. Another source of the discrepancy is that dif- 
ferent studies assume different wavelength dependences 
for the emissivity. In general when a shallower emissivity 
function is chosen, this produces warmer temperatures. 

In this paper we combine far-infrared ISO photometry 
with submillimeter photometry from the Submillimetre 
Common-User Bolometer Array (SCUBA) (Holland et 
al. 1999) at the James Clerk Maxwell Telescope (JCMT) 
to sample the full spectrum of cool dust emission in nor- 
mal spiral galaxies. We will first present the dust temper- 



atures we infer from 60, 100, and 180 /<m ISO measure- 
ments for 71 spiral galaxies in the ISO Atlas of Bright 
Spiral Galaxies (Bendo et al. 2002) (henceforth referred 
to as Paper 1). We then examine dust temperatures ob- 
tained by combining the far-infrared data with submil- 
limeter data for 8 of the galaxies in the sample. Finally, 
we calculate dust masses and examine gas to dust mass 
ratios for 7 of these galaxies. 

2. Data 

2.1. The Sample 

The galaxies in this sample are a subset of a com- 
plete, magnitude-limited set of galaxies selected from 
the Revised Shapley-Ames (RSA) Catalog (Sandage & 
Tammann 1987). The sample comprised galaxies with 
Hubble types between SO and Sd and with magnitudes 
Bt = 12 or brighter; galaxies in the Virgo Cluster were 
excluded. A randomly selected, subset of these galaxies 
was observed by ISO based on target visibility. This pro- 
duced a total of 77 galaxies that are representative of the 
range of Hubble types in the RSA Catalog, 71 of which 
were observed successfully with ISOPHOT. Detailed in- 
formation on the sample is presented in Paper 1. 

2.2. Observations and Data Analysis 

2.2.1. Far-Infrared Photometry 

The 60, 100, and 180 (im photometric data were taken 
with ISOPHOT (Lemke et al. 1996) and processed with 
PIA 8.0 (Gabriel et al. 1997) with additional processing 
to correct for PSF effects, as discussed in Paper 1. Note 
that, for the analysis without the submillimctcr data, wc 
assume that the 180 /xm flux within the 180 arcsec aper- 
tures effectively comes from the same sources as the 60 
and 100 /im fluxes within a 135 arcsec aperture. To inves- 
tigate how much extra far-infrared flux may come from 
the region outside the 135 arcsec aperture but within the 
180 arcsec aperture, we investigated the distribution of 
infrared flux in maps of IRAS data made by the HIRES 
software package provided by the Infrared Processing and 
Analysis Center. We find that the median ratio of flux 
within 135 arcsec to flux within 180 arcsec is 0.85 at 
60 /xm and 0.80 at 100 ^m for the galaxies in this sam- 
ple. According to the HIRES beam profiles, point sources 
would have ratios of 0.95 at 60 ^m and 0.85 at 100 fim. 
The small difference between the ratio of fluxes for the 
sample galaxies and the ratio expected for point sources 
suggests that most of the emission between the 135 and 
180 arcsec apertures in the HIRES maps is mostly from 
the point spread function and that emission in those re- 
gions is relatively small in most galaxies. We will assume 
that the distribution of the 180 [im emission will be sim- 
ilar to the 100 /im emission, so therefore, we will assert 
that most of the 180 ^m emission for the galaxies in this 
sample also comes from within 135 arcsec. When we also 
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have submillimeter data, however, we can make assump- 
tions about the spatial distribution of the emission and 
perform a deconvolution analysis to estimate the 180 ^im 
fluxes within 45 and 135 arcsec apertures. 

2.2.2. Submillimeter Photometry 

We observed NGC 4088 and NGC 4826 at 850 /mi 
using SCUBA at the JCMT. The data were taken in jig- 
gle map mode, in which the telescope's secondary mirror 
shifts in a 64-point pattern to create a completely sam- 
pled map for the inner 135 arcsec of the target. For sky 
subtraction, we nodded 2 arcmin off-source parallel to 
the optical disks' minor axes for each galaxy. Integra- 
tion times were 64 min for NGC 4088 and 60 min for 
NGC 4826. We performed skydips regularly to check the 
opacity during the observations. At 850 /um, the opacity 
typically ranged between 0.3 and 0.6 for all the observa- 
tions. (The high sky opacity made the simultaneous 450 
//m observations useless.) For calibration, we observed 
Mars using similar methods. 

In addition to the above data, we have used data 
from the JCMT SCUBA archives for 6 galaxies: NGC 
4414, NGC 4631, NGC 5236, NGC 5713, NGC 5792, and 
NGC 5907. All of these data were taken in SCUBA's 
jiggle map mode, with cither Uranus, OH 231.8, or 
IRC +10216 used for calibration. 

The data were processed with the SCUBA User Re- 
duction Facility (Jenness & Lightfoot 1998) following the 
typical processing steps. First, the data were flat-fielded 
and corrected for airmass and sky opacity. Bad bolome- 
ters were flagged and masked. An additional background 
selection was made by measuring and subtracting surface 
brightnesses in selected off-target bolometers using the 
"remsky" procedure. Finally, the data were combined 
together and rbinned to produce final maps. Images of 
planets and other calibration sources were processed us- 
ing the same methods. Signals were measured within 
45 and 135 arcsec apertures, which are comparable to 
the apertures of ISOPHOT. For calibration, we measured 
the signals from the standard sources in similarly-sized 
apertures, and these then provided conversion factors for 
converting signals to flux densities for the targets. We as- 
sume that the uncertainty in both the new measurements 
and the archival measurements comes primarily from the 
uncertainty in the calibration, which is approximately 
10 %. Other sources of uncertainty, mainly uncertainties 
from background noise, appear to be less significant. Ta- 
ble 1 lists the resulting flux densities as well as the optical 
diameter parameter from the Third Reference Catalogue 
of Bright Galaxies (de Vaucouleurs et al. 1991) and an 
indication of whether the data came from original obser- 
vations or the JCMT SCUBA archives. 



3. Spectral Energy Distributions and Dust Tem- 
peratures 

3.1. Far-Infrared Spectral Energy Distributions 

First, we fit blackbody functions with two different 
emissivity functions to the data. The emissivity func- 
tions we selected scale as A -1 and A~ 2 , functions that 
are used by many studies cited in Section 1 and that are 
recommended by Hildebrand (1983). 

Table 2 lists the dust temperatures determined from 
the ISO far-infrared data, and Figure 1 shows histograms 
of these temperatures. We estimate an error of ~ 2 K in 
the dust temperatures due to uncertainties in the pho- 
tometric calibration of ISOPHOT. This error was deter- 
mined by fitting blackbodics to data where the 60 /tin 
flux was increased 20 % and the 180 /tm flux was de- 
creased 20 % and vice-versa. The temperatures of the 
fits are, for the the A -1 emissitivity, a mean of 31 K with 
a range of 24 - 42 K, and, for the A~ 2 emissitivity, a mean 
of 25 K with a range of 21 - 32 K. The functions with the 
A -1 emissitivities generally fit the data better according 
to the x 2 criterion. 

The dust temperatures vary little in relation to mor- 
phology and AGN activity. The data in Table 3 do show 
that some SO and elliptical galaxies tend to have warmer 
dust temperatures, but all the spiral and irregular galax- 
ies have indistinguishable dust temperatures. 

These results should be compared with those of Stickel 
et al. (2000), who used IRAS 60 and 100 /mi data supple- 
mented by 170 /xm photometry from the ISO Serendip- 
ity Survey for 115 galaxies. They assumed an emis- 
sivity function scaling as A~ 2 and found a temperature 
distribution centered at 20 K with a range of 15-25K. 
These results are systematically colder by 5-6 K and 
they conclude that "for the first time this indicates that 
a cold dust component with < 20 K is present in 
a large fraction of all spiral galaxies." In contrast, us- 
ing the same emissivity function we find no galaxy with 
Tp < 21 K and no support for such a cold dust com- 
ponent in most spiral galaxies. One possible explanation 
could be that Stickel et al. (2000) use IRAS Faint Source 
Catalog (FSC) data at 60 and 100 /im, whereas we have 
used ISO photometry at 60 and 100 /im and were able to 
make approximate corrections for the point spread func- 
tion. As we demonstrated in Paper 1, the ISO fluxes used 
in this research, particularly the 60 /im fluxes, are sys- 
tematically higher than those in the IRAS FSC Catalog, 
which would naturally lead to the temperatures found 
using the ISO photometry to be higher than those found 
using corresponding IRAS FSC data. 

However, some qualifications must be made to inter- 
pretations based on data confined to this relatively small 
wavelength range, even though these data from ISO in- 
clude the longer wavelength 180 /im measurement which 
IRAS did not have. First, if dust components 10-15 K ex- 
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ist, their contribution would not be evident using the far- 
infrared data alone. Second, emission at 60 /im may in- 
clude emission from warmer dust components that cause 
the spectral energy distribution to flatten out (see Lau- 
reijs et al. (2000) for an example based on observational 
data of galactic cirrus, Popcscu ct al. (2000) and Dunne 
& Ealcs (2001) for examples of blackbody fits to ob- 
servations of nearby galaxies, or Desert, Boulanger, & 
Puget (1990) and Dale et al. (2001) for detailed modeling 
of dust emission at far- infrared wavelengths). Thus, fit- 
ting single blackbody functions through the 60, 100, and 
180 /im data points is possible but may not accurately de- 
scribe the physics of the dust emission. The addition of 
longer wavelength data may provide evidence for more 
than one dust component. Therefore, in the next sec- 
tion, we examine the blackbody fits obtained when the 
far-infrared data are combined with submillimeter data. 

3.2. Far-Infrared to Submillimeter Spectral En- 
ergy Distributions 

We have submillimeter and far-infrared data for a sub- 
set of 8 of the ISO galaxy sample, and blackbody func- 
tions with A -1 and emissivities, with (3 variable, 
have been fitted to these data. We refer to the resulting 
dust temperatures measured within 135 arcsec apertures 
as global dust temperatures. We also determined the 
temperatures of dust within the inner 45 arcsec of the 
galaxies, and we refer to these as the nuclear dust tem- 
peratures. (We note our designation "nuclear" is very 
qualitative; at the distances of these galaxies 45 arcsec 
subtends > 1 kpc, a much larger region than would nor- 
mally be considered the galaxy nucleus. Starburst re- 
gions in galaxies are typically 250 - 500 pc.) The dust 
temperatures measured within a 135 arcsec aperture but 
outside a 45 arcsec aperture we refer to as disk temper- 
atures. These distinctions permit investigation of the 
differences between dust temperatures in the inner and 
outer regions of these galaxies. 

Since the diffraction limit of the ISO telescope at 
180 /im is ~ 145 arcsec, to do this analysis for the 180 /im 
data requires estimating the fluxes from both central un- 
resolved and extended components. We use a method 
similar to that used by Radovich et al. (1999). (We also 
used this method for the 60 and 100 /im data in Paper 1.) 
We calculate the ratio of fluxes for the point and extended 
sources at both 100 /im and 850 /im, since we have the 
required 45 arcsec resolution at both these wavelengths. 
We then interpolate to estimate the ratio of the 180 /im 
fluxes from point and extended components. According 
to Laurcijs (1999) 14 , 

61.2 % of the flux from a point source centered in the 
C200 array (the 180 /im measurement) will fall on the 
array. Therefore the flux within the 180 arcsec aperture 



can be expressed as 

/ m (180") - 0.612/ c + (180") 2 .f e 



(1) 



where f m is the measured flux, f c is the flux from the 
central source, and f e is the flux per arcsec 2 from a dif- 
fuse, uniform, extended source. Substituting the ratio 
- — jf- — into the above equation allows for solving for f c 

and f e : 
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(3) 



These central and extended components can then be used 
to calculate fluxes within 45 and 135 arcsec apertures: 



/ t (45")-/ c +(45") 2 /e 



/ t (135") = / c + (135") 2 /, 



1 http://www.iso. vilspa.esa.es/users/expl_lib/PHT/c200fpsf02.ps.gz 



(4) 
(5) 

The resulting 180 /im flux desnities are listed in Table 4. 
We can now use these 180 /im flux estimates for the dust 
temperature analysis outlined above. 

Table 5 lists the resulting dust temperatures derived 
using the A -1 emissivities, and Table 6 lists the dust 
temperatures for A _/3 emissivities. In addition, Figures 2 
- 9 show the spectral energy distributions and the fits 
of blackbody functions with \~P emissivities. We also 
tried fitting single blackbodies with A -2 emissivities to 
the data, but except for NGC 5236 (which will be dis- 
cussed below) the fits were very poor. We estimated 
an error of ^ 1 K in these dust temperatures using the 
same method described above for the far-infrared spec- 
tral energy distributions. In summary, the fits yield dust 
temperatures of ~ 35 K for models with A -1 emissivities 
and ~ 30 K for models with A _/3 emissivities. Typical 
fits to the emissivity coefficient were between 0.9 and 1.4. 

NGC 5236 is exceptional in that the emissivity func- 
tion is closer to A~ 2 . NGC 5236 is a galaxy with unusual 
supernova activity. Shocks from the strong supernova 
activity associated with the galaxy's nucleus could be re- 
sponsible for breaking up larger dust grains into smaller 
grains, thus making the wavelength dependence of the 
emissivity function steeper. The disk of NGC 4414, but 
not the nucleus, also exhibits a steep emissivity. Star 
formation has been inhibited in the nucleus by the lack 
of molecular gas present (Sakamoto 1996), but sites of 
strong star formation are evident in the disk (c/. Fig. 25 
in Paper 1), and the large dust grains may have suffered 
similar fragmentation as in NGC 5236. 

In general, the "nuclear" dust temperatures are sys- 
tematically warmer than the disk temperatures in all 
cases. Because stellar densities are higher in the nuclei of 
these galaxies, and because the star formation rate may 
be higher, the interstellar radiation field is stronger. This 
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enhances both the dust emission and the dust tempera- 
tures, producing the observed differences in dust temper- 
atures between the nuclear and disk regions. 

Since some other studies have suggested there are 
two temperature components, we have also investigated 
the possibility that the dust emission may be composed 
of two blackbody components. This obviously requires 
steeper wavelength-dependence for the emissivities. We 
therefore chose emissivity functions proportional to A~ 2 
for both temperature components. Table 7 shows the 
best fitting temperatures, with two temperatures derived 
for all galaxies except NGC 5236. (For NGC 5236, the 
analysis demonstrates that a second blackbody compo- 
nent would make only a minor contribution to the far- 
infrared and submillimeter spectral energy distribution, 
so only one component was fit to the data.) In most 
cases, the two components have temperatures of ~ 15 K 
and ~ 30 K. For galaxies with no 450 /an data, the so- 
lution is an even-determined problem with an exact fit 
that should be interpreted with caution. Moreover, the 
colder dust temperature is dependent mainly on the sub- 
millimeter measurement, so errors in the submillimeter 
measurement directly affect the temperature of the cold- 
est component. 

Having used submillimeter data to determine dust 
temperatures for these 8 galaxies, we can now compare 
the temperatures found for these same galaxies used the 
ISO far-infrared data alone. Comparing the dust tem- 
peratures in Table 2 with the "Global" temperatures in 
Table 5 and Table 7 shows that the addition of the sub- 
millimeter data systematically increases the derived dust 
temperatures ~ 2-3 K. The more significant difference is 
that, by adding the submillimeter data point(s) one can 
search for a very cold dust component of temperature 
- 6-15 K. 

In summary, for these 8 galaxies with both ISO far- 
infrared and submillimeter measurements, we find dust 
temperatures in almost all cases that can be described 
as a ~ 35 K blackbody with a A -1 emissivity, a ~ 30 K 
blackbody with a variable emissivity coefficient ranging 
from 0.9 to 1.4, or two blackbodies with temperatures of 
~ 15 K and ~ 30 K and A~ 2 emissivities. Clearly the far- 
infrared-to-submillimeter spectral energy distributions of 
galaxies can be represented with several different func- 
tions and correspondingly larger or smaller dust temper- 
atures. This is essentially a duplication of the variety of 
results found in the literature for the past two decades. 
The key question is, which of these is most likely to be 
correct? We return to this question in the Discussion. 

4. Gas to Dust Mass Ratios 

The accurate measurement of dust temperatures is not 
only necessary for understanding dust emission in galax- 
ies. It is also required for measuring the gas to dust 
ratios in spiral galaxies. As discussed by Devereux & 



Young (1990), the temperatures and dust masses mea- 
sured with IRAS data, when combined with molecular 
gas data, produced gas to dust mass ratios of approxi- 
mately 1000. This is a factor of 10 larger than the value 
for the Milky Way and it therefore raises a serious prob- 
lem. This problem has also appeared in most other stud- 
ies in which cold dust components have been found. The 
ISO Scrcndipidy Survey by Stickel et al. (2000) is a no- 
table exception; the typical gas-to-dust mass ratios found 
are near that of the Milky Way, but the ratios range from 
50 to 1000. If these results are reliable they raise the 
question whether this basic property of the interstellar 
medium really varies so dramatically among galaxies. 

Young et al. (1995) list CO flux measurements for the 
8 galaxies for which we have far-infrared and submillime- 
ter flux measurements. Using their data and their con- 
version of CO intensities to molecular gas surface densi- 
ties, we calculated molecular gas densities. We then used 
the 850 /im flux densities and the dust temperatures and 
the associated emissivity assumptions used to the fit the 
SEDs to calculate dust masses. Working at 850 ^m is 
preferred because the determination of dust masses at 
that wavelength is relatively insensitive to uncertainties 
in dust temperatures. With these masses, we then cal- 
culated gas to dust ratios that may be compared to the 
Milky Way ratio. 

4.1. Dust Masses and Gas to Dust Mass Ratios 

We calculated dust masses using 

r F(v)D 2 lr 4a . 
M Dust = [ ) J }[——}p 6 

derived in Hildcbrand (1983). is the dust mass, F(z/) 
is the measured flux density, D is the distance, B(v,T) 
is the flux density from a blackbody of temperature T, 
a is the approximate average size of an interstellar dust 
grain, Q(v) is the emissivity, and p is the density of a dust 
grain. Since these are nearby spiral galaxies, we assume 
that these galaxies contain dust with properties similar to 
the dust in the Milky Way. Therefore, we have also used 
the asssumptions in Hildebrand (1983) for calculating 
dust masses. The size and density of the dust grains are 
assumed to be a = 0.1 ^m and p = 3 g cm~ 3 . For the 
emissivities, we use the expression 

Q(A,/?) = Q a(^V (7) 

from Klaas et al. (2001), where Q = 40 cm -1 and (3 is the 
emissivity coefficient. This equation gives Q(850^m,l) = 
1.2 x 10 L4 and Q(850,um,2) = 3.5 x 10~ 5 . 

Table 8 lists the global dust masses calculated from 
the data in Tables 1, 5, 6, and 7. Dust masses calcu- 
lated with steeper emissivity laws are always higher than 
those calculated with shallower emissivity laws for two 
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reasons. First, the fits to the spectral energy distribution 
using the shallower emissivity function produce higher 
temperatures than those fits using the steeper emissivity 
functions. These higher temperatures lead to higher val- 
ues of B(i/,T), which then result in lower calculated dust 
masses. The second is the difference in the dust grain 
emissivity values at 850 /an. The emissivity is higher for 
shallower emissivity functions, which leads to lower cal- 
culated dust masses. Essentially, a higher temperature 
or a higher emissivity requires less dust to produce the 
measured submillimeter flux. 

Table 8 also presents global molecular gas masses cal- 
culated from the data in Young et al. (1995) and the re- 
sulting global gas to dust mass ratios. (We did examine 
gas to dust ratios for separate nuclear and disk regions, 
but the comparison showed no clear distinction between 
nuclear and disk ratios.) 

4.2. Discussion 

Inspection of Table 8 shows that the gas to dust ra- 
tios ultimately depend on the emissivity function used to 
fit the data. Perhaps even more importantly, when the 
emissivity coefficient is fixed (i.e. proportional to A -1 or 
A~ 2 ), the range of gas to dust ratios for these 8 galaxies 
varies by factors of 20 - 30. With the A -1 emissivities, the 
ratios range from 130 to 2700, and with the A~ 2 emissiv- 
ities and two temperature components, the ratios range 
from 13 to 480. However, when the emissivity coefficient 
is fit to the data, the gas to dust ratios vary from 150 to 
580, only a factor <~ 4. 

We may check the plausibility of these calculated gas 
to dust ratios using the Milky Way ratio as a benchmark. 
The molecular gas to dust ratio in the Milky Way has 
been determined to be 115 using COBE data (Sodroski 
et al. 1994). The median gas to dust ratio for the A -1 
emissivity fits is approximately a factor of 6 higher than 
the ratio for the Milky Way, so the single blackbody with 
the A -1 emissivity is unlikely to accurately represent the 
physical conditions of the dust. The median ratio from 
the A -2 emissivity fit is within a factor of 3 of the Milky 
Way value, so these fits are more plausible. The me- 
dian ratio from the emissivity fits is also within a 
factor of 3, so these fits are equally plausible. Dale & 
Helou (2002), however, have demonstrated empirically 
that fitting one or two blackbodies to far-infrared and 
submillimeter data may lead to underestimates of the 
dust mass by as much as a factor of 10 when compared 
to more realistic models that fit a series of blackbodies 
to the data. If we accept that this effect is present, the 
A -1 and A _/3 emissivity fits would produce gas to dust 
ratios near the Milky Way's but the A -2 emissivity ratios 
would be much too low to be plausible. 

However, the range in ratios provides an additional 
criterion for comparing the three fits. The dust masses 
computed using both the A -1 emissivity and the two- 



component model with very cold dust give gas-to-dust 
ratios which vary from galaxy to galaxy by factors of 20 
and 37 respectively. In contrast, the range in the ratios 
for the A _/3 fits is less than a factor of 4. Since these 
are all nearby spiral galaxies, we would expect the basic 
properties of the interstellar medium to be similar. If 
the best description of the dust is the one that gives the 
lowest dispersion in the gas-to-dust ratios, and gas/dust 
ratios which are similar to that for the Milky Way, then 
the \~P emissivity fits are clearly preferred. The fits 
with the fixed emissivities give implausible variations in 
the gas-to-dust ratios, while the A _/3 fits give a more 
acceptable range. 

As a result, using both closeness to the Milky Way gas- 
to-dust ratio and the small dispersion of the ratios as an 
indication of the plausibility for the gas-to-dust ratios, it 
is evident that the fits using the A _/3 emissivity, i.e. the 
right-hand column of Table 8, provide results that are 
most likely to be correct for the galaxies in this sample. 

We are now able to address the larger question of the 
lack of consensus in the dust temperatures and even the 
number of different temperature components in normal 
galaxies which was described in the Introduction. It is 
clear from the analysis above that one can fit spectral 
energy distributions from 60 - 850 /im equally well by a) 
a single temperature component with emissivity propor- 
tional to A -1 , b) a single temperature component with 
emissivity proportional to A - * 3 , or c) two temperature 
components, one of which is "cold," with a temperature 
~ 15 K, and the other with temperature ~ 30 K. This 
range of results is just what one finds in reviewing the 
literature on this subject over the past 20 years. We 
suggest it is the assumptions made about the emissivity 
function in different studies which have largely produced 
the wide range of conclusions about the dust tempera- 
tures and dust masses in galaxies and the lack of consen- 
sus which has plagued this subject. We have shown that 
that the assumptions most commonly used lead to im- 
plausible gas-to-dust ratios. The distinguishing feature 
of this study is that we have tested the physical plausi- 
bility of various possible emissivity functions rather than 
adhering rigidly to an assumed emissivity function and 
calculating dust temperatures and other parameters. 

So, having shown that the emissivity fits to the 
spectral energy distributions imply more plausible dust 
masses and gas-to-dust ratios than the other choices, 
what is the physical interpretation of this result? First 
we note that in this model we use only one tempera- 
ture component to approximate the dust emission when 
there is almost certainly a range of dust temperatures 
both within star formation regions and in the diffuse in- 
terstellar medium. Does this make sense? A priori, a 
two-component model might be more reasonable since 
a warmer component from star formation regions and a 
cooler component from dust heated by the diffuse inter- 
stellar radiation field are expected. Since such models 
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appear to result in dust masses that are too large and 
gas-to-dust ratios that are implausibly small, this sug- 
gests that the colder component in such models is a less 
significant contributor than the two-component fits im- 
ply. Since only a single temperature gives plausible dust 
masses, we suggest that if a second cold component to 
the dust emission from the diffuse interstellar medium is 
present, it is relatively insignificant in terms of both its 
contributions to the mass and to the SED. To distinguish 
its contribution apparently requires better spectral cov- 
erage and photometric precision in the submillimeter and 
millimeter regions than has been available heretofore. 

We must also consider the physical interpretation of a 
variable emissivity law due to the varying j3 we find for 
this galaxy sample. The emissivity functions should be 
related to the average physical properties of the dust, 
such as the composition and grain size distribution, 
and these are expected to be relatively constant among 
nearby spiral galaxies. Is it not inconsistent to argue 
that the emissivity functions can vary when we also ar- 
gue that the gas-to-dust ratios should be similar among 
spiral galaxies? We suggest two plausible interepreta- 
tions for this result. First, it may be that the grain 
properties do vary among these galaxies. Dale & Hclou 
(2002) and references therein present data that demon- 
strate how j3 may depend on environment. However, 
they argue that j3 should decrease as the temperature 
of dust increases, which is the opposite of the trend we 
have found. Instead, we suggest that galaxies which have 
undergone significant recent star formation activity and 
concomitant supernovae may have a somewhat altered 
dust grain population. Supernova shock waves can break 
down silicate grains with carbon mantles into individual 
silicate and graphite grains. Silicate grains have differ- 
ent emissivities than the graphite grains, so the average 
emissivity would change. Furthermore, the destruction of 
very large grains would the make the emissivity function 
steeper. Breaking down the dust with supernova shock 
waves also suggests some dust grain destruction (Dyson 
& Williams 1997), which is exacly what is seen in these 
results for galaxies with steep emissivities. An alterna- 
tive interpretation might be that the varying emissivity 
function is simply a second fitting parameter, in addition 
to a temperature, and this combination better accounts 
for the range of dust temperatures and other dust prop- 
erties that are found in spiral galaxies. 

5. Conclusions 

We have shown that, for the 71 galaxies in the ISO At- 
las sample, the 60 - 180 /im spectral energy distributions 
have mean temperatures of 31 K for the A -1 cmissitivi- 
ties and 25 K for the A -2 emissitivities. The blackbodies 
with A -1 emissitivities generally fit the data better than 
the ones with A -2 emissivities. Analysis of the dust tem- 
peratures as a function of morphological type reveals that 



spirals from Sa-Scd have remarkably uniform dust tem- 
peratures which are virtually identical within the errors. 
Dust temperatures are enhanced above the average for 
the spirals only in E/SO galaxies and in galaxies with 
especially strong star formation activity. Data based on 
SEDs over this spectral range do not support the exis- 
tence of a substantial cold component at a temperature 
< 20 K in a large fraction of spiral galaxies, in contrast 
to some earlier studies. 

For the eight galaxies in the sample with the spec- 
tral energy distributions measured from 60 - 850 ^m, we 
showed that the 60 - 850 /im spectral energy distribu- 
tions could be fit equally well by: a) a single tempera- 
ture component with emissivity proportional to A -1 , b) 
a single temperature component with emissivity propor- 
tional to A - ^, or c) two temperature components, one of 
which is "cold," with a temperature ~ 15 K. We suggest 
that the enormous variation in reported dust tempera- 
tures and dust masses that characterize the literature on 
spiral galaxies is largely a reflection of the variety of as- 
sumptions for the emissivity functions employed. 

To choose among these three assumptions we used the 
gas-to-dust ratio as a criterion. Only the temperatures 
derived using the single temperature component with 
emissivity proportional to produced gas-to-dust ra- 
tios within a rather narrow range and one that was close 
to the gas-to-dust ratio for the Milky Way. We suggest 
that these results are most likely to represent the actual 
dust properties of normal spiral galaxies. In particular, 
we find the wide range in inferred gas-to-dust ratios de- 
rived with the assumption of a "cold" dust component 
less physically plausible than those ratios derived using 
a single temperature component with temperatures de- 
rived using A _/ % with (3 calculated for each galaxy. 
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Fig. 1.- 
data. 



Histograms of the temperatures with (a) A 1 and (b) A 2 cmissivities, as determined from the 60 - 180 jim 
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Fig. 2. — The (a & b) global, (c & d) nuclear, and (e & f) disk spectral energy distributions for NGC 4631. The measured 
data are represented by squares, and the errors for the data are less than the size of the squares. Parts a, c, and e show 
the best fitting blackbodies with A _/3 emissivities (with j3 variable). Parts b, d, and f show the best fitting two blackbody 
components with A~ 2 emissivities, with the dotted lines represent the individual blackbody components. This format also 
applies to Figures 3-9, although only the global spectral energy distributions are shown in those cases. 
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Fig. 3. — The global spectral energy distribution for NGC 4088. 




Fig. 4. — The global spectral energy distribution for NGC 4414. 




Fig. 5. — The global spectral energy distribution for NGC 4826. 
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Fig. 6. — The global spectral energy distribution for NGC 5236. 




Fig. 7. — The global spectral energy distribution for NGC 5713. 




Fig. 8. — The global spectral energy distribution for NGC 5792. 
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Fig. 9. — The global spectral energy distribution for NGC 5907. 
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Table 1 

submillimeter flux density measurements 



Galaxy 450 /um Flux 850 /um Flux log log D 2 5 a Data 

Density (Jy) Density (Jy) 





45 arcsec 


135 arcsec 


45 arcsec 


135 arcsec 




Source 


NGC 4088 






0.26 ±0.03 


0.98 ±0.10 


1.76 


Observation 


NGC 4414 


2.8 ±0.3 


7.3 ±0.7 


0.36 ±0.04 


0.84 ±0.08 


1.56 


Archive 


NGC 4631 


4.8 ±0.5 


18. ±0.2 


0.54 ±0.05 


1.89 ±0.19 


2.19 


Archive 


NGC 4826 






0.63 ±0.06 


1.01 ±0.10 


2.00 


Observation 


NGC 5236 


5.5 ±0.5 


9.7 ±0.9 


0.82 ±0.08 


1.36 ±0.14 


2.11 


Archive 


NGC 5713 






0.26 ±0.03 


0.43 ± 0.04 


1.44 


Archive 


NGC 5792 






0.18 ±0.02 


0.33 ±0.03 


1.84 


Archive 


NGC 5907 






0.43 ± 0.04 


1.26 ±0.13 


2.10 


Archive 



a The optical diameter parameter defined in de Vaucouleurs et al. (1991). 
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Table 2 

Temperatures from Far-Infrared Spectral Energy Distributions 



Galaxy Temperature (K) a 

A -1 emissivity A~ 2 cmissivity 



NGC 55 


26 


22 


NGC 289 


31 


25 


NGC 1512 


31 


26 


NGC 3359 


32 


26 


NGC 3556 


31 


26 


NGC 3898 


30 


25 


NGC 4062 


28 


23 


NGC 4088 


31 


25 


NGC 4096 


29 


24 


NGC 4100 


33 


26 


NGC 4136 


30 


25 


NGC 4157 


29 


24 


NGC 4203 


32 


26 


NGC 4236 


32 


26 


IV T / 1 /"i J C\ A A 

NGC 4244 


27 


23 


IV T / 1 /"i A C\r-t A 

NGC 4274 


29 


24 


NGC 4314 


34 


28 


NGC 4395 


26 


22 


TV T/"1 /~1 A A 1 A 

NGC 4414 


31 


25 


NGC 4448 


30 


25 


NGC 4559 


29 


24 


NGC 4605 


32 


26 


NGC 4618 


30 


25 


NGC 4631 


34 


27 


NGC 4710 


33 


27 


NGC 4725 


24 


21 


NGC 4826 


33 


27 


NGC 4984 


39 


30 


NGC 5033 


30 


25 


NGC 5054 


31 


26 


at/"( /^n r /~i r r 

NGC 5055 


28 


24 


NGC 5087 


38 


30 


NGC 5101 


28 


24 


NGC 5102 


29 


24 


NGC 5112 


29 


24 


NGC 5170 


28 


24 


NGC 5204 


33 


27 


NGC 5236 


36 


29 


NGC 5247 


29 


24 


NGC 5300 


29 


24 


NGC 5334 


28 


23 


NGC 5364 


26 


22 


NGC 5371 


28 


23 


NGC 5457 


27 


23 


NGC 5474 


32 


26 


NGC 5556 


33 


27 
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Table 2 — Continued 



Galaxy Temperature (K) a 

A -1 emissivity A~ 2 emissivity 



NGC 5566 


29 


24 


NGC 5584 


32 


26 


NGC 5585 


33 


27 


NGC 5669 


29 


24 


NGC 5676 


31 


26 


NGC 5701 


31 


25 


NGC 5713 


37 


29 


NGC 5746 


28 


23 


NGC 5792 


33 


27 


NGC 5838 


42 


32 


NGC 5846 


37 


29 


NGC 5850 


30 


25 


NGC 5866 


32 


26 


NGC 5907 


26 


22 


NGC 5985 


28 


23 


NGC 6015 


29 


24 


NGC 6215 


35 


28 


NGC 6217 


35 


28 


NGC 6221 


34 


27 


NGC 6300 


31 


26 


NGC 6340 


38 


30 


NGC 6503 


30 


25 


NGC 6643 


31 


26 


NGC 6744 


25 


21 


NGC 6753 


32 


26 



a Errors on the temperature measurements are 
2 K. 
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Table 3 

Dependence of Far-Infrared Dust Temperatures on Morphological Type 



Morphological 
Type 


Number 


Mean Temperatures (K) 
with A -1 Emissivity 


Mean Temperatures (K) 
with A -2 Emissivity 


All 


71 


30.9 ± 0.4 


25.3 ± 0.3 


E - SO/a 


11 


34.5 ± 1.4 


27.5 ± 0.8 


Sa - Sab 


8 


30.0 ± 1.1 


25.0 ± 0.8 


Sb - Sbc 


19 


30.3 ± 0.7 


24.9 ± 0.5 


Sc - Scd 


25 


30.2 ± 0.5 


24.9 ± 0.3 


Sd - Irr 


8 


30.9 ± 1.1 


25.4 ± 0.8 
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Table 4 

Results of Deconvolution Analysis on 180 /zm Data 



Measured 180 /im Deconvolved 180 /xm Flux Densities (Jy) 
Galaxy Flux Densities in Components in Set Apertures 

(180") (Jy) Central Extended (135") (45") (135") 



NGC 4088 


46.2 


8.4 


23. 


11. 


31. 


NGC 4414 


54.5 


22. 


23. 


25. 


45. 


NGC 4631 


93.9 


19. 


46. 


24. 


65. 


NGC 4826 


51.8 


18. 


34. 


34. 


50. 


NGC 5236 


192. 


90. 


77. 


99. 


170. 


NGC 5713 


22.3 


16. 


7. 


17. 


23. 


NGC 5792 


15.9 


8.7 


6.0 


9.3 


15. 


NGC 5907 


44.1 


9.1 


22. 


12. 


31. 
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Table 5 

Temperatures with A -1 Emissivity 



Galaxy Global Nuclear Disk 





Temp. (K) a 


Reduced x 2 


Temp. (K) a 


Reduced x 2 


Temp. (K) a 


Reduced x 2 


NGC 4088 


33. 


0.014 


35. 


0.44 


32. 


0.17 


NGC 4414 


34. 


3.0 


35. 


0.95 


33. 


12. 


NGC 4631 


35. 


3.4 


37. 


3.5 


34. 


3.7 


NGC 4826 


35. 


2.0 


36. 


2.2 


32. 


2.5 


NGC 5236 


45. 


10. 


48. 


7.0 


41. 


18. 


NGC 5713 


39. 


1.6 


41. 


2.9 


34. 


0.64 


NGC 5792 


35. 


0.79 


37. 


1.4 


30. 


1.26 


NGC 5907 


27. 


0.33 


28. 


0.46 


27. 


0.26 



a Errors on the temperature measurements are 2 K. 
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Table 6 

Temperatures with \~ f) Emissivity 



Galaxy 




Global 




Nuclear 




Disk 






Temp. 





Reduced 


Temp. 





Reduced 


Temp. 


P 


Reduced 




(K) a 




2 

X 


(K) a 




2 

X 


(K) a 




2 

X 


NGC 4088 


34. 


1.0 


0.0055 


33. 


1.1 


0.11 


34. 


0.9 


0.0030 


NGC 4414 


29. 


1.4 


0.52 


31. 


1.3 


0.66 


26. 


1.7 


5.4 


NGC 4631 


35. 


1.0 


5.1 


34. 


1.2 


3.8 


35. 


0.9 


5.9 


NGC 4826 


31. 


1.3 


0.011 


32. 


1.3 


0.11 


28. 


1.4 


0.20 


NGC 5236 


30. 


1.9 


0.19 


33. 


1.7 


0.80 


26. 


2.1 


0.69 


NGC 5713 


34. 


1.3 


0.022 


35. 


1.4 


0.15 


31. 


1.2 


0.32 


NGC 5792 


32. 


1.2 


0.28 


34. 


1.2 


0.82 


27. 


1.3 


0.17 


NGC 5907 


28. 


0.9 


0.45 


29. 


0.9 


0.83 


28. 


0.9 


0.28 



a Errors on the temperature measurements are 2 K. 
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Table 7 

Temperatures with A~ 2 Emissivity 



Galaxy Global Nuclear Disk 

Temp. (K) a Reduced x 2 Temp. (K) a Reduced y 2 Temp. (K) a Reduced y 2 



NGC 4088 


13 


& 


29. 




14. 


& 30. 




12 


& 28. 




NGC 4414 


15 


& 


28. 


0.031 


7. 


& 27. 


0.70 


18 


& 28. 


6.1 


NGC 4631 


12 


& 


30. 


2.5 


14 


& 31. 


1.9 


12 


& 30. 


2.9 


NGC 4826 


14 


& 


28. 




15 


& 30. 




5. 


& 25. 




NGC 5236 




29. 




0.29 




30. 


1.6 




27. 


0.78 


NGC 5713 


14 


& 


31. 




16 


& 33. 




6. 


& 27. 




NGC 5792 


15 


& 


30. 




17 


& 34. 




9. 


& 24. 




NGC 5907 


13 


& 


26. 




14 


& 27. 




13 


& 25. 





a Errors on the temperature measurements are 2 K. 
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Table 8 

Global Dust Masses and Gas to Dust Ratios 



Galaxy 


Dust Masses (M 


3 ) 




Molecular Gas 




Gas / Dust Ratios 






A- 1 


A- 2 




8 


Masses (M ) 


A- 1 


A- 2 


A _/3 




cmissivity 


emissivity 


emissivity 




emissivity emissivity emissivity 


NGC 4088 


4.8 x 10 6 


5.7 x 10 7 


4.6 x 


10 6 


1.5 x 10 9 


310. 


26. 


330. 


NGC 4414 


1.3 x 10 6 


1.2 x 10 7 


2.6 x 


10 6 


1.5 x 10 9 


1200. 


130. 


580. 


NGC 4631 


1.4 x 10 6 


2.1 x 10 7 


1.4 x 


10 6 


4.4 x 10 8 


310. 


21. 


310. 


NGC 4826 


2.7 x 10 5 


2.7 x 10 6 


4.5 x 


10 5 


1.8 x 10 8 


670. 


67. 


400. 


NGC 5236 


3.5 x 10 5 


2.0 x 10 6 


1.7 x 


10 6 


9.5 x 10 8 


2700. 


480. 


560. 


NGC 5713 


5.5 x 10 6 


6.5 x 10 7 


9.3 x 


10 6 


3.0 x 10 9 


550. 


46. 


320. 


NGC 5907 


6.2 x 10 6 


5.9 x 10 7 


5.2 x 


10 6 


7.8 x 10 s 


130. 


13. 


150. 
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